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Abstract 



We present results of searches for i?-meson decays to charmless final states 
that include a K*, p, uj, or (p meson accompanied by a second meson. Using the 
entire data sample of 9.7 x 10^ BB pairs collected with the CLEO II and CLEO 
II. V detectors, we observe a signal for the decay B^ — s- ujtt~^, and measure a 
branching fraction of B{B+ ujtt~^) = (ll.Sllg ± 1.5) x 10~^. We also see 
evidence for the decay B^ — > luK^, and set limits for the decays B^ p^n^ and 
B^ — > K*^TT^. In addition to these new results, we also summarize previous 
CLEO results on related channels. All quoted results are preliminary. 
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I. INTRODUCTION 



The study of charmless hadronic decays of B mesons plays a key role in understanding 
the phenomenon of CP violation within the Standard Model Large asymmetries are 
predicted for some exclusive final states, and although the relatively small branching fractions 
of O{10~^) currently limit the experimental reach for such studies, the sensitivity of the 
CLEO II detector allows us to measure branching fractions for some decay modes. B factory 
experiments that just started or are about to begin operation should expand this already 
rich field of investigation. 

Theoretical predictions typically make use of effective Hamiltonians, often with factor- 
ization assumptions P-[TT|. The strong interaction between particles in the final state com- 
plicates these predictions. However, experimental measurements can be used to verify the 
validity of the assumptions made, to tune the parameters of the theory in order to make fur- 
ther predictions, and to understand the relative importance of the various decay amplitudes 
that contribute to a particular decay. Recently, it has been suggested |T2|-p^ that published 
experimental results on charmless hadronic B decays indicate that cos 7 < 0, where 7 is 
one of the angles of the unitarity triangle. This somewhat disagrees with current fits to the 
information most sensitive to CKM matrix elements [0. Again, more experimental results 
can help clarify the situation. 

In this paper, we present preliminary results of searches for 5-meson decays to exclusive 
two-body final states |]T3[ that include a K*, p, u, or meson, and another low- mass me- 
son. We concentrate on new results, although previously reported results are included for 
completeness. 



II. DATA SAMPLE AND INITIAL EVENT SELECTION 



The data were collected with the CLEO II fT^ and CLEO II. V detectors at the 



Cornell Electron Storage Ring (CESR). For most of the new results the data sample includes 
all data collected prior to the de- commissioning of CLEO II. V, in preparation for a significant 
upgrade. The total integrated luminosity is 9.13 fb~^ for the reaction e'^e~ — > T(45') — > BB, 
which corresponds to 9.7 x 10^ BB pairs. This is between 40% and a factor of three more 



statistics than for previously published results [[19i ^. In addition, we re-analyzed the CLEO 



II data set with improved calibration constants and track-fitting procedure, allowing us 
to extend our geometric acceptance and track quality requirements. This has led to an 
overall increase in reconstruction efficiency of 10-20 % compared to the previously published 
analyses. For studies of background from continuum processes, we also collected 4.35 fb~^ 
of data at a center-of-mass energy below the threshold for BB production. 

The final states of the decays under study are reconstructed by combining detected 
photons and charged pions and kaons. The resonances in the final state are identified via the 
decay modes p — > tttt, K* — > Ktt, uj — >■ 7r"'"7r~7r°, and (f) — > K^K~ . The detector elements 
most important for the analyses presented here are the tracking system, which consists of 
several individual concentric detectors operating inside a 1.5 T superconducting solenoid, 
and the high- resolution electromagnetic calorimeter, made of 7800 CsI(Tl) crystals. For 
CLEO II, the tracking system consists of a 6-layer straw tube chamber, a 10-layer precision 
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drift chamber, and a 51-layer main drift chamber. The main drift chamber also provides 
a measurement of the specific ionization loss, used for particle identification. For CLEO 
II. V the 6-layer straw tube chamber was replaced by a 3-layer double-sided-silicon vertex 
detector, and the gas in the main drift chamber was changed from an argon-ethane to a 
helium-propane mixture. 

Reconstructed charged tracks are required to pass quality cuts based on their track fit 
residuals and impact parameter in both the r-0 and r-z planes, and on the number of main 
drift chamber measurements. Each event must have a total of at least four good charged 
tracks. The specific ionization [dE/dx) measured in the drift layers is used to distinguish 
kaons from pions. Expressed as the number of standard deviations from the expected value, 
5*4 (i = TT or K), it is required to satisfy 15*^1 < 3.0. Electrons are rejected based on 
dE/dx and the ratio of the measured track momentum and the associated shower energy in 
the calorimeter. Muons are rejected by requiring that charged tracks penetrate fewer than 
seven interaction lengths of material. Pairs of charged tracks used to reconstruct K^s (via 
TT+TT^) are required to have a common vertex displaced from the primary interaction 
point. The invariant mass of the two charged pions is required to be within two standard 
deviations (10 MeV) of the nominal mass. Furthermore, the momentum vector, 
obtained with a mass-constrained kinematic fit of the charged pions' momenta, is required 
to point back to the beam spot. 

Photons are defined as isolated showers, not matched to any charged tracks, with a lateral 
shape consistent with that of photons, and with a measured energy of at least 30 (50) MeV 
in the calorimeter region | cos 6*1 < 0.71(> 0.71), where 9 is the polar angle. Pairs of photons 
are used to reconstruct vr^s. The momentum of the pair is obtained with a kinematic fit of 
the photons' momenta with the 7r° mass constrained to its nominal value. The resolution of 
the invariant mass of the two photons depends on the momentum of the 7r° and is between 
5 and 10 MeV/c^. We require the reconstructed mass to be within 3(7 on the low side of the 
nominal 7r° mass, and 2cr on the high side. 



III. ANALYSIS TECHNIQUE 

The primary means of identification of B meson candidates is through their measured 
mass and energy. The quantity /S.E is defined as AE" = Ei + E2 — Ei,, where Ei and E2 
are the energies of the two mesons in the final state, and Eb is the beam energy. The beam- 
constrained mass of the candidate is defined as M = El — |pP, where p is the measured 
momentum of the candidate. We use the beam energy instead of the measured energy of the 
B candidate to improve the mass resolution by about one order of magnitude. 

The large background from continuum quark-ant i- quark (gg) production can be reduced 
with event shape cuts. Because B mesons are produced almost at rest, the decay products 
of the BB pair tend to be isotropically distributed, while particles from qq production have 
a more jet-like distribution. The cosine of the angle 9$ between the sphericity axis |^ 



of the charged particles and photons forming the candidate B and the sphericity axis of 
the remainder of the event should have a fiat distribution for B mesons and be strongly 
peaked at ±1.0 for continuum background. We require | cos 6*5! < 0.8. For final states 
containing an or meson we use the thrust axis instead of the sphericity, and we 
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require |cos^t| < 0.9. We also form a Fisher discriminant (JF) with the momentum 
scalar sum of charged particles and photons in nine cones of increasing polar angle around 
the sphericity axis of the candidate, the angle of the sphericity axis of the candidate with 
respect to the beam axis, and R2 = H2/H0, the ratio of the second and zeroth Fox- Wolfram 
moments [^. For analyses with an or in the final state, the thrust axis replaces the 
sphericity axis, and the angle between the candidate's momentum vector and the beam axis 
replaces i?2- 

The specific final states investigated are identified via the reconstructed invariant masses 
of the B daughter resonances. For vector-pseudoscalar final states, further separation of sig- 
nal events from combinatoric background is obtained through the use of the defined angular 
helicity state of the vector meson in the final state. The observable 7i is the cosine of the 
angle between the direction of the B meson and the vector meson daughter decay direction 
(normal to the decay plane for the u), both in the vector meson's rest frame. 

Signal event yields for each mode are obtained with unbinned, multi-variable maximum 
likelihood (ML) fits. We also perform event counting analyses that apply tight constraints 
on all variables described above. Results for the latter are consistent with the ones presented 
below. For each input event, the likelihood (£j) is defined as 

m 
i=l 

where Vi are the probabilities for each of the m hypotheses of the fit, and n^, the free 
parameters of the fit, are the number of events in the overall sample for each hypothesis. 
The Vi are the product of the probability distribution functions (PDFs) for each of the fit's 
input variables. For input events, the overall likelihood is then 

1=1 

where the first term takes into account the Poisson fluctuations in the number of events. In 
all cases, the fit includes hypotheses for signal decay modes and the dominant continuum 
background. For a few channels indicated below, we also include a hypothesis for background 
from other B decay modes. For all others, we verified that this component is negligible. 

The variables used in the fit are AE, M, JF, resonance masses, and Ti as appropriate. 
For pairs of final states differentiated only by the identity (charged pion or kaon) of one of 
the two mesons, we also use the dE/dx measurement, S*,, for the high- momentum track and 
fit for both modes simultaneously. Correlations between input variables were investigated 
and found to be negligible. For each decay mode investigated, the signal PDFs for the input 
variables are determined with fits to high-statistics Monte Carlo event samples generated 



with a GEANT based simulation of the CLEO detector response. The parameters of 
the background PDFs are determined with similar fits to the sum of off-resonance data and 
a sideband region of on-resonance data in Ai? and M. For M, the sideband is defined by 
5.2 < M < 5.3 GeV/c^. For AE the sideband varies depending on the decay mode; details 
are given in the sections below. Sideband regions for each of the other input variables are 
also included in the likelihood fit sample. 
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IV. RESULTS 



Table | gives all the measurement results. Details specific to various final states are given 
in separate sections below. Results for decay modes with a (j) meson in the final state and 
of the types p^h^ and K*^h^ , where is a charged pion or kaon, are based on a data 
sample of 5.8 x 10^ BB pairs. Results for final states p~h^, K*^h^, and p^K^ are based 
on a sample of 7.0 x 10^ BB pairs. Shown in the table are the signal event yield, the 
efficiency, the product of the efficiency and the relevant branching fractions of particles in 
the final state, the statistical significance of the observed yield in standard deviations, the 
branching fraction central value with statistical and systematic error, and the corresponding 
90% confidence level upper limit. The one standard deviation statistical error on the central 
value is determined by finding the values where the quantity = ~21n(£/£max) changes 
by one unit, where £max is the point of maximum likelihood. 

Systematic errors are separated into two major components. The first is systematic er- 
rors in the PDFs, which are determined with variations of the PDF parameters within their 
uncertainty, taking into account correlations between parameters. The second component 
is systematic errors associated with event selection and efficiency factors. These are de- 
termined with studies of independent data samples. For branching fraction central values, 
the systematic error is the quadrature sum of the two components. For upper limits, the 
likelihood function is integrated to find the yield value that corresponds to 90% of the total 
area. This value is then increased by its systematic error, and the efficiency is reduced by 
one standard deviation of its systematic error when calculating the final upper limit. 

Table H shows the final results for each decay mode investigated. The third column 
of the table indicates whether the result is from this work or from an earlier analysis. For 
observations the final result is reported as a branching fraction central value, while for modes 
where the yield is not sufficiently significant we quote the 90% confidence level upper limit. 
Also in the table are previously published theoretical estimates. 

A. Final states including an u meson. 

The results for decay modes including an uj meson in the final state were obtained with 
the full data sample (9.7 x 10^ BB pairs). The final selection prior to the likelihood fit 
requires \AE\ < 200 MeV. In Table |, the final state uh'^ represents the sum of the ujK~^ 
and uj7r~^ states = K'^ or tt^). For uK* and up final states, the resonance mass sidebands 
are 830 < Mk* < 950 MeV/c^ and 600 < Mp < 950 MeV/c^, respectively, and cross-feeds 
between those decay modes are ignored when performing the fit. For K*^ — >■ K^tt", there are 
two possible assignments for the K^. The choice is made based on the dE/dx information 
of the tracks. For p° ti^ti^ the two charged tracks are simply assumed to be pions. 

The maximum likelihood fit for ujK*^ and up^ includes a hypothesis for background from 
generic B decays. PDFs for this hypothesis are obtained from a sample of Monte Carlo 
generated events that corresponds to approximately 1.7 times the size of the on-resonance 
data sample. The yield for up'^ is entirely consistent with cross-feed from uK*^. The quoted 
branching fraction central value has been adjusted to take this into account. For the final 
states u!K*~^ and cup"*", the 7r° from K*~^ or p+ decay defines the daughter direction. Since the 



7 



TABLE I. Measurement results. Columns list the final states with resonance decay modes as 
subscripts, event yield from the fit, reconstruction efficiency e, total efficiency including secondary 
branching fractions Bs, statistical significance (a), branching fraction central value B, and the 
corresponding 90% confidence level upper limit. For central values the first error is statistical and 
the second systematic. 



Final state 


Yield(events) 


e{%) 


eBs{%) 


Signif. 


B{10-^) 


90% UL(10 


ujK+ 

L0h+ 

ujK*+ 


28.5+^-^ 

i.5i?:i 

7.9li? 

10.811:1 
3.7+1-0 

i.olfl 
7.011:^ 


29 
22 
29 
24 
29 

7.1 
18 

6.8 
14 


26 
19 
26 

7.4 
26 

6.3 
16 

2.0 

8.3 


6.2 
0.6 
2.1 
3.9 
7.3 
2.8 
0.9 
0.3 
2.3 


ii.3i^:i^ ± 1.5 
o.sij:^ ±0.5 
3.2i^:^ ±0.8 
lo.oil:^ ± 1-5 
i4.3i|:^ ± 2.1 

181^^ ± 6 

n n+5-7 +2.9 
"-'•"-'-0.0 -0.0 

Btf ± 6 

9.iil:J ± 1.9 


17 
5.8 
8.0 
21 
21 
47 
11 
52 
19 


p^K+ 
p-K+ 
p^K^ 


26.11IJ 

ZO.O_j g 

o 4+5:2 

14.8+^-? 
8.3llig 

Q 0+4-9 


30 
12 
34 
28 
11 


30 

12 
34 
28 
11 
10 


5.2 
5.6 

2.7 


1511 ± 4 
35+|i ± 5 


5.1 
22 
25 
27 




12.31^-? 

o.il^:? 
o.olg:J 


37 


18 
24 
18 






27 
4.2 
12 


K*+n- 


io.8lti 

5.7lti 




7 
4.1 


5.2 
2.5 
5.9 


23l? ± 3 
9n+i5 +3 

99+8 +4 
^^-6 -5 




KkXoK- 
K*+K- 


o.olH 




7 
4.1 


0.0 
0.0 
0.0 




8 

17 
6 



07rO 

<pK+ 2.4lf0 
<PK^ 4.3l|2 



54 27 0.0 

35 17 0.0 

53 26 1.3 

41 7.0 2.6 



4.0 
5.4 

1.6li:^±0.2 5.9 
10.7l^:^ ± 1.1 28 
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TABLE II. Final results and expectations from theoretical models. 



Decay mode 




Source 


Theory S (10 ^) 


7~l-l- 


4- 


11.3^29 =t 


This work 


0.6 - 11 


B^' - 


n 


< 5.8 


This work 


0.01 - 12 


- 




< 8.0 


This work 


0.2 - 13 


B^ - 




< 21 


This work 


0.02 - 10 


B^ - 




14.3^3:2 ± 2.1 


This work 




B^ - 




< 47 


This work 


7-28 


B^ - 


n 


< 11 


This work 


0.005 - 0.4 


B^ - 




< 52 


This work 


0.9 - 15 






< 19 


This work 


0.3 - 12 


B+ - 




15l| ± 4 


m 


0.4-8 


B^- 




35l}J ± 5 


m 


26-52 


B^- 




< 5.1 


This work 


0.9 - 2.3 


B+ - 


^p'>K+ 


< 22 


m 


0.1 - 1.7 




p-K+ 


< 25 


m 


0.2 - 2.5 


B^- 


-^p'^K'^ 


< 27 


m 


0.04-1.7 


B+ - 




< 27 




4-12 


B^- 




99+8 +4 
^^-6 -5 


m 


1.2 - 19 


B^- 




< 4.2 


This work 


1.1 - 5 


B^- 


^ K*+K- 


< 6 


m 




B+ - 


K*^K+ 


< 12 


m 


0.2 - 1 


B+ - 


-> (j)7r~^ 


< 4.0 


m 


0.001 - 0.4 


B^- 




< 5.4 


m 


0.0004 - 0.2 


B+ - 


(/)K+ 


< 5.9 


m 


0.3 - 18 


B^- 




< 28 


m 


0.3 - 18 



References 
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distribution of Ti. is not known for these vector-vector final states we assume the worst case 
(Ti.^) when computing the efficiency, and we require 7i < to reduce the large combinatoric 
background from soft vr^'s. This also results in much reduced cross-feed between the two 
decay modes. 

As can be seen in Table |, we observe a clear signal for the decay un^. Figure |I| 

shows the contours in Icr intervals of the likelihood function. Figures ^ and ^ show projections 
of the data on fit variables after tight selection cuts were applied on the variables other than 
the one plotted. The solid curves are not a fit to the data showed in each plot, but rather 
an overlay of the fit function scaled to take into account the additional cuts applied. 




0.0 0.5 1.0 1.5 2.0 

B(B+^aDK^) (10"^) 



FIG. 1. Likelihood function contours for — > utt^ and — > ujK^ . Systematic errors are 
not included in the contours. 



B. The decays B^ p^ir^ and B^ K*^-k^. 

The results for the decay modes p°7r° and B^ — ^ K*^7t^ are based on a sample of 

8.3 X 10^ BB pairs. There is significant cross-feed between the two decay modes, due in 
large part to the presence of a high momentum 7r° in the final state, which results in poorer 
AE resolution. As a first step both decay modes are fit simultaneously. In this case, event 
selection prior to the maximum likelihood fit requires the invariant mass of the resonance to 
be within 0.3 < Mh+h- < 1-0 GeV/c^, as well as -0.3 < AE < 0.2 GeV. Both quantities 
are computed assuming that the two charged tracks are pions. The continuum background 
hypothesis in the ML fit is separated into four components, depending on the identity of 
the charged tracks (vr^ or K^). For input to the ML fit, AE, the resonance mass, and Ti 
are calculated assuming the two charged tracks are pions for the p°7r° signal and continuum 
background hypotheses, and assuming one charged track to be a kaon and the other a pion 
for the K*^iT^ hypothesis. We also use the dE/dx information for the two charged tracks. 
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□ 





5.200 



COTI^ fit 

Background 




5.225 5.250 5.275 

B mass (GeV/c^) 



5.300 



> 



o 



a 

CD 
> 



12 



10 



-0.20 



fit 

Background 



-0.10 



0.00 



0.10 



AE (GeV/cO 



0.20 



FIG. 2. Projection onto the reeonstrueted B mass (left) and AE (right) for — > ujtt^. The 
sohd hne shows the result of the likelihood fit, scaled to take into account the cuts applied to 
variables not shown. The dashed line shows the background component. 



(OtC fit 

Background 




0.73 0.755 0.78 0.805 

7i"^7i'7i° mass (GeV/c^) 



0.83 



□ COJt* 

am* fit 

Background 



0.00 




0.25 



0.50 

Fisher 



0.75 



1.00 



FIG. 3. Projection onto the reconstructed u mass (left) and J- (right) for B'^ — >■ o^tt"*". The 
solid line shows the result of the likelihood fit, scaled to take into account the cuts applied to 
variables not shown. The dashed line shows the background component. 
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The results of this combined fit estabhshed that there is no significant signal for the decay 
mode p°7r°. For this mode, the combined fit results are reported in Tables | and ||. 

A second ML fit was then performed for K*^7r'^ only. The preliminary selection 

required the dE/dx information for both the charged pion and charged kaon candidates in 
the final state to be within 2a of the expected value. We also required |A£'| < 0.2 GeV, and 
that the resonance mass be within 150 MeV/c^ of the known K*^ mass. The ML fit looked 
for only two hypotheses, signal K*^ti^ and continuum background. In this fit, no dE/dx 
information was used as input. No significant yield was observed, as indicated in Table |. 



V. DISCUSSION 

The observation of the decay B^ —>■ uj-k^ yields a branching fraction in complete agree- 
ment with the previously reported measurement for B^ p^n^. For the corresponding 
neutral decay modes B^ — > un^ and B^ —>■ p^vr^, fairly stringent upper limits are set, as 
could be expected because of their color suppression relative to the charged decay modes. 
For the latter decay mode, the upper limit from this work is another piece of information 
for studies of decays of the type B pvr, to go along with the previous observations of 
B+ p%+ and B° ^ p-n+ @. 



We also see evidence for the decay mode B^ ujK^. In this case the charged decay 
mode B^ — >■ ujK^ is expected to have a similar branching fraction. For this latter mode, 
the additional data and re-analysis of old data did not support the previously reported 



observation |]T9l. However, the central value for B^ ujK^ is only about l.So" away from 
the central value for B^ —>■ uK^ . There is no significant evidence for decays of the type 
B pK, although the results are consistent with the ones for B —>■ uK decays. 

The analysis presented here finds no significant yield for the decay mode B^ — > K*^tt^, 
in agreement with theoretical expectations which point to a branching fraction smaller than 



the previously observed B^ — > K*^'k ||2^. There is no evidence for decay modes of the type 



B K*K, as well as for (f)K and 07r final states. 



VI. CONCLUSIONS 

We have observed the decay mode B~^ — >■ uj7t~^, and measure a preliminary branching 
fraction of B{B^ — > uott^) = (11.3^|^±1.5)xl0-^ in agreement with the previously reported 
observation of B^ p^n^. We also see evidence for the decay B^ uoK^. New upper 
limits are set for other final states including an uj meson, and for the decays B^ p^n^ 
and K*'^tt^. In combination with previously reported measurements and upper limits, 

these results are generally in agreement with predictions based on factorization models in 
charmless hadronic B decays. 
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